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ABSTRACT 

A cold dark matter paradigm predicts that galaxy clusters follow an universal mass density 
profile and fit a well defined mass-concentration relation, with lensing clusters being pref- 
erentially triaxial haloes elongated along the line of sight. Oddly, recent strong and weak 
lensing analyses of clusters with a large Einstein radius suggested those haloes to be highly 
over-concentrated. Here, we investigate what intrinsic shape and orientation an halo should 
have to account for both theoretical predictions and observations. We considered a sample 
of 10 strong lensing clusters. We first measured their elongation assuming a given mass- 
concentration relation. Then, for each cluster we found the intrinsic shape and orientation 
which are compatible with the inferred elongation and the measured projected ellipticity. We 
distinguished two groups. The first one (nearly one half) seems to be composed of outliers 
of the mass-concentration relation, which they would fit only if they were characterised by 
a filamentary structure extremely elongated along the line of sight, that is not plausible con- 
sidering standard scenarios of structure formations. The second sample supports expectations 
of iV-body simulations which prefer mildly triaxial lensing clusters with a strong orientation 
bias. 
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1 INTRODUCTION 

Clusters of galaxies, the most recent bound structures to form 
in a hierarchical cold dark matter model with a cosmological 
constant (ACDM), offer important clues to the assembly pro- 
cess of structure in the universe. A^-body simulations are success- 
ful in fitting large-scale structure measurements and are able to 
mak e detailed theoretical predictions on dark matter halo proper- 
ties dNavarro et alj 19971 : IBullock etalj200lt iDiemand et alj2004l ; 



iDuffv et aL 20081) butsome disagreement with observations still 
persists. One possible conflict between ACDM and measurements 
is the detec tion of extremely large Einstein radii in massive lens- 
ing cluster jBroadhurst & Barkanj|2008l:ISadeh & Rephaelil l 2008l : 
IOguri&Blandforj|2009l : IZitrinetal.ll2009l) . The Einstein radius 
mirrors the mass contained in the inner regions and its measurement 
is quite model indepe ndent. Even if an universa l Navarro-Freank- 
White density profile dNavarro et al.1 199611 1991 NFW) reproduces 
many characteristics of massive lenses, such haloes should be over- 
concentrated to fit the data. 

The concentration parameter measures the halo central den- 
sity, which depends on the assembly history and thereby on the 
time of formation. The halo concentration is then expected to be 
related to its viri al mass, with the co ncentration decreasing grad- 
ually with mass IBull ock etal]|200lh . Concentrations of massive 
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galaxy clusters are then a crucial probe of the mean density of the 
universe at relatively late epochs. State-of-the art models of cos- 
mic structure formation suggest that galaxy cluster concentrations 
decrease gradually with virial mass. However, cluster observations 
have yet to firmly confirm this correlation. 

On the observational side, the situation at present is 
still unclear due to t he p lurality of methods employed 
dComerford & Nataraiar] 120071) , The observed concentration- 
mass relation for galaxy clusters has a slope consistent with 
theoretical prediction from si mulations, though the nor malization 
factor seems to be higher (Comerford & Nataraiar] 1 20071) . A 
critical point is that concentrations measured in massive lensing 
clusters appear to be systematically larger than X-ray concen- 
trations dComerford & NataraiarJ |2007|) . A simil ar, though less 
prono unced, effect is also found in simulations dHennawi et al.l 
120071) . which show that massi ve lensing clusters a r e usua lly elon- 
gated along the line of sight. Oguri & Blandford (2009) showed 
that the larger the Einstein radius, the larger the over-concentration 
problem, with clusters looking more massive and concentrated due 
to the orientation bias. 

The over-concentration bias seems to be much la rger in obser- 
vations than in simulations. Broadhu rst et alj d2008l) inferred sig- 
nificantly high concentrations for four nearly relaxed high-mass 
clusters. Such a trend has been recently exacerbated with the anal- 
ysis of the largest known Einstein rad ius in MACS J07 17.5+3745 
dZitrin et alj|2009h . lOguri et alj d2009h found that the data from a 
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sample of ten clusters with strong and weak lensing features were 
highly inconsistent with the predicted concentration parameters, 
even including a 50% enh ancement to account for the lensing bias 
(Oguri&Blandford 2009). On the other hand, Ok abe et al, (2009) 
found that the correlation in the c — M relation, as measured in 
a sample of 19 clusters with significant weak lensing signal that 
were well fitted by a NFW profile, was marginally compatible with 
predictions for both slope and normalization. 

Different definitions of parameters for spherically averaged 
profile s can play a role when comparing observations to predic- 
tions dBroadhurst & Barkanal l2008h . Triaxiality issues were ad- 
dressed by ICorless et al. I d2009hT who derived weak lensing con- 
straints on three strong lensing clusters without assuming a spheri- 
cal halo model. The large errors that accompany triaxial parameter 
estimates can make observations compatible, even if marginally, 
with theoretical predictions. Investigations in the weak lensing 
regime demonstrated that neglecting halo triaxiality can lead to 
over- and under-estimates of up to 50 % and a factor of 2 i n halo 
mass and concentration, respectively dCorless & Kindl2007l) . An 
analysis of AC 1 14 using only strong lensing data and accounting 
for triaxiality also supported that projecti on effects play a m ajor 
role in the estimate of the concentration (Sereno et al. 200^). Fi- 
nally, analyses of stacked weak lensing cluste rs of lesser mass does 
not exhibit the high conc entration problem dJohnston et alj|2007l ; 
iMandelbaum et al1l2008l). in agreement with theoretical findings 
dOguri & Blandfordll2009l) . 

Several effects can play a r ole: over-concentrated clusters have 
a larger lensing cross section dHennawi et ail [2007 ) ; strong lens- 
ing clusters pr eferentially sample the high-m ass end of the cluster 
mass function (Come rford & Nataraianl2007h ; while extreme cases 
of triaxiality are rare, such halos can be much more efficient lense s 
than their more spherical counterparts dOguri & Blandfordl l2009); 
the strongest lenses in the universe are expected to be a highly 
biased population preferentially orientated alon g the line of sight 
dHennawi et alj200l;lo"g~uri & Blandfo rd 2009); estimates of lens- 
ing concentrations can be also inflated due to substructures close 
to the line of sight dPuchwein & Hilbert 2009). On the other hand, 
contamination of weak lensing catalogues c an lead to underesti- 
mate the concentration dLimousin et al .120071) . 

In order to check the ACDM paradigm is then crucial to 
account for all possible biases when comparing theoretical re- 
lations with lensing observa tions. Such approach was taken in 
iBroadhurst & Barkanal d2008l) . who derived the probability distri- 
bution of Einstein radii from concentration distributions found in 
iV-body simulations. Also after considering that lensing clusters 
are intrinsically over-concentrated and that the inherent triaxiality 
of CDM haloes along with the presence of substructure enhances 
the projected mass in some orientations, they fo und that theoretical 
predic tions are excluded at a 4a significance. Sadeh & Rephaeli 
d2008l) reached a similar conclusion. They implied the Einstein ra- 
dius distribution from the probability distribution of cluster forma- 
tion times and from a formation redshift-concentration scaling de- 
rived from iV-body simulations. However due to various inherent 
uncertainties, the statistical range of the predicted distribution may 
be significantly wider than commonly acknowledged. 

Here, we compare measurements with theoretical predictions 
from semi-analytical investigations and iV-body simulations avoid- 
ing some possible biases connected to spherical averaging. The 
paper is as follows. In Section [2] we review the predictions from 
either iV-body simulations or semi-analytical investigations. Sec- 
tion [3] discusses how projected quantities are related to intrinsic 
parameters for an ellipsoidal cluster. In Sec. [4] we develop our in- 



version method which under some given a-priori hypotheses allows 
to infer intrinsic mass, concentration and elongation of a lensing 
cluster; the method is then applied to a sample of ten strong lensing 
clusters. In Sec.|5]we compare the observed distributions of elonga- 
tion along the line of sight and ellipticity in the plane of the sky to 
different theoretical predictions. Section [6] exploits the previously 
inferred geometrical parameters to predict the intrinsic axial ratios 
and the orientation of the clusters in the sample. Finally, Sec.|7]is 
devoted to a summary and to some final considerations. 

Throughout the paper, we assume a flat ACDM cosmology 
with density parameters Qm =0.3, Qa — 0.7 and an Hubble con- 
stant Ho = 100ft km s _1 Mpc _1 , h — 0.7. We quote uncertainties 
at the 68.3% confidence level. 



2 THEORETICAL PREDICTIONS 

High resolution A-body simulations have shown that the density 
profiles of dark matter halos are successf ully described as NFW 
density profiles dNavarro et al ]fl996llT997h . whose 3D distribution 
follows 



Pnpw 



Ps 



(r/r s )(l + r/r s )2 : 



(1) 



where p s is the characteristic density and r s is the characteristic 
length scale. A-body simulations showed as well that haloes are 
aspherical and that such profiles can be accu rately described b y 
concentric triaxial ellipsoids with aligned axes (Jing & Sutol2002h . 
A NFW equivalent profile whose density is constant on a family 
of similar, concentric, coaxial ellipsoids is obtained by replacing 
the spherical radius r with an ellipsoidal radial variable £ in the 
intrinsic orthogonal framework centred on the cluster barycentre 
and whose coordinates, Xi,j n t, are aligned with its principal axes, 



c 2 



i=l 



(2) 



where ei are the intrinsic axial ratios. Without loss of generality, 
we can fix ei ^ e% ^ e$ = 1. In the following, we will also use 

the inverse ratios, < qi = 1/e; 1. 

According to recent TY-body simulations ( Neto et al J [20071 : 
iMaccio et al]|2008l ; iGao et al.ll2008l: iDuffv et alj|2008l) . the depen- 
dence of dark matter halo concentration c on halo mass M and 
redshift z can be adequately described by a power law 



c = A(M/M pivot ) B (l + z) c 



(3) 



Since several assumptions were used by competing groups, re- 
sults can be somewhat different, in particular as far as the over- 
all normalization is concerned. Several values for the linear am- 
plitude of mass fluctuations as were considered. The higher as, 
the earlier the formatio n epoch for haloes of a given mass. Here, 
we follow IDuffv et al] d2008h , who used the cosmological pa- 
rameters from WMAP5 (as = 0.796) and found {A,B,C} = 
{5.71 ± 0.12, -0.084 ± 0.006, -0.47 ± 0.04} for a pivotal mass 
Mpi vo t = 2 x 10 12 Mqj/h in the redshift range — 2 for their full 
sample of clusters. 

By separately studying the distribution of NFW profile param- 
eters both for the general halo population and for the lensing pop- 
ulation (i.e. ha l oes w eighted by their strong lensing cross-section) 
iHennawi et al.| j2007) showed that the distribution of 3D concentra- 
tions of the lens population is the same as that of the general halo 
population except for a shift upwards by a factor of ~ 17%. In the 
following, we will then also consider an enhanced c — M relation 
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for lensing clusters, with A ~ 6.68. Note that such increased value 

of A could be also seen as due to a larger value of <jg. 

jV-body simulations prefer mildly triaxial halos. Ijing & Sutd 

J2002h investigated the probability distribution of intrinsic axial ra- 
tios and proposed an universal approximating formula for the dis- 
tribution of minor to major axis ratios, 



200 



C200 



P(qi) oc exp 



(gi - q»/r qi ) 2 
2al 



where g M = 0.54, a. 

r n = (Af t , ir /M.)°- orn « W 



0.113 and 

0.7 



(4) 



(5) 



with M* the characteristic nonlinear mass at redshift z. The condi- 
tional probability for 92, the ratio of the intermediate to the major 
axis-length, goes as 



P(qi/q2\qi) = 



2(1 



2 qi /q 2 - 1 



1 



(6) 



for 91/52 ^min = max[gi, 0.5], whereas is null otherwise. The 
lensing population has nearl y the same triaxiali ly distribution as 
the total cluster population jHennawi et al J 1 20071) . This could be 
explained as the result of two counter-balancing effects. Whereas 
both triaxiality and concentration increase the lensing cross section, 
the shape of a dark halo is correlated with its concentration, with 
more concentrated clusters being more spherical. 

For comparison we will also consider a flat distribution for the 
axial ratios, such that 



1 



for the full range < qi ^ 1 and 

P(g 2 |gi) = (1-gi)" 1 



(7) 



(8) 



for Q2 ^ qi and zero otherwise. The resulting probability for q? is 
then P(q2) = ln(l — 92)" . Such a flat distribution allows also 
for very triaxial clusters (qi < q 2 <C 1), which are preferentially 

excluded by TV-body simulations. 

Finall y, semi-analitycal dOguri & Blandford l2009h and nu- 
merical taennawi etal] 12007) investigations showed a large ten- 
dency for lensing clusters to be aligned with the line of sight. De- 
noting the angle between the major axis and the line o f sight as 8, 
such condition can be expressed as dCorless et ail2009ri 



P(cos 9) oc exp 



(cos 61 - l) 2 



(9) 



with ae = 0.115. For comparison, we will also consider a popula- 
tion of clusters randomly oriented, i.e. 

P(cos0) = l (10) 
forO < cos (9 «S 1. 

3 PROJECTION OF TRIAXIAL HALOES 

Dealing with ellipsoidal halos, we need g eneralized definit i ons fo r 
the intrinsic NFW parameters. We follow [Corless & Kind d2007l) . 
who defined a triaxial virial radius r^oo such that the mean density 
contained within an ellipsoid of semi-major axis T2oo is A = 200 
times the critical density at the halo redshift; the corresponding 
concentration is C200 = ?"200 / r s ■ Then, the characteristic overden- 
sity is expressed in terms of C200 as for a spherical profile, 



3 ln(l + C200) - C2oo/(l + C200) 



(11) 



The virial mass, M200, is the mass within the ellipsoid of semi- 
major axis r2oo, M200 = (8007r/3)gig2r2ooPc r . Such defined C200 
and A/200 have small deviations with respect to the parameters 
computed fitting spherically averaged density profiles, as done in 
TV-body simulations. The only caveat is that the spherical mass ob- 
tained in simulations is significantly less than the ellipsoidal M200 
for extreme axial ratios dCorless & Kina2007l) . However, since the 
dependence of the concentration on the mass is quite weak, see 
Eq. <[3j, this will have negligible effects on our analysis. 

Three rotation angles relate the intrinsic to the observer's coor- 
dinate system, i.e. the three Euler's angles, 8, ip and ip. After align- 
ment of the observer's coordinate system with the direction con- 
necting the observer to the cluster centre, the line of sight has polar 
angles {8,p — 7r/2} in the intrinsic system. With a third rotation, 
ip, we can properly align the coordinate axes in the plane of the sky. 
If not stated otherwise we will line up such axes with the axes of 
the projected ellipses. 

When viewed from an arbitrary direction, quantities co nstant 
on sim ilar ellipsoids project themselves on similar ellipses JStarkl 
In general, the projected map F20 on the plane of the sky 
and the intrinsic spheroid al volume density P3D are related by 
dStarklll977r . lSerendl2007h . 



F 2 T>(£;lp,Pi) = -j= 



C 



(12) 



where £ is the elliptical radius in the plane of the sky, l s is the 
typical length scale of the 3D density, Ip is its projection on the 
plane of the sky, pi are the other parameters describing the intrinsic 
density profile (slope, ...) and / is a function of the cluster shape 
and orientation, 



2 2 2 2 2 

/ = ex sin 8 sin cp + e 2 sin 1 



2 2 
cos ip + cos ( 



(13) 



the subscript P denotes measurable projected quantities. 

Let us see in some details how the parameters describing the 
projected map depend on the intrinsic shape and orientation of the 
3D distribution. The axial ratio of the major to the minor a xis of the 
obser ved projected isophotes, ep(^ 1), can be written as (Binggeli 
1 19801) . 



ep 



1 3 + 1 + y/(j-Q 2 +4fc 2 

' j + i- V(i-0 2 +4fc 2 ' 

where j, k and I are defined as 



2 2-2 

3 = e^sin 1 



.2 2 , 
+ ei cos ( 



cos 2 p + el cos 2 1 



• 2 

sin ip, 



k — (ei — e 2 ) sin tp cos p cos 8, 
I 



2 • 2 
e-i sin ip ■ 



2 2 
e 2 cos p>. 



(14) 



(15) 
(16) 
(17) 



In the following we will also use the ellipticity e = 1 — 1/ep . 

The observed scale length Ip is the projectio n on the plane of 
the sky of the cluster intrinsic length dStarkl 19771) . 



Z D — li 



(— ) 

V eie2 / 



1/2 



1/4 



(18) 



(19) 



Equation {18) can be rewritten as 

Is _ Ip 

where the parameter eA quanti fies the elon gation of the triaxial 
ellipsoid along the line of sight |Sereno 20q3), 
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6A = 



Veie2/ 



1/2 



r 



3/4 



(20) 



The quantity lp /e& represents the half-size (along the line of sight) 
of the ellipsoid as seen from above, i.e. perpendicularly to the line 
of sight. If eA < 1, then the cluster is more elongated along the 
line of sight than wide in the plane of the sky. The smaller the eA 
parameter, the larger the elongation. In the following, we will use 
as an elongation parameter also a geometrical factor 



/geo 



/ Nl/2 1/2 

(eiea) ' e t ' 



/3 /4 



eA 



(21) 



Summarising, the surface density can be expressed in terms of 
projected quantities as 

F-2D = —f2o(^,ep,-4>;lp;p i ,...), (22) 
eA 

where /2D has the same functional form as for a spherically sym- 
metric halo. Iln order to write Eq. d22t in its actual form, we ex- 
ploited that the integral in £ in Eq. H2\ is proportional to the in- 
trinsic scale length l s . The dependence on the elongation eA is 
decoupled from the dependence on the apparent ellipticity and in- 
clination. The other parameters characterising the 3D profile only 
account for the radial dependence of the projected density. Then, 
when we deproject a surface density, the normalization of the vol- 
ume density can be known only apart from a geometrical factor. 
Note that in our notation, the elliptical radius is written as a func- 
tion of the coordinates in the plane of the sky as 



£ 2 = {x\ + epxl)(l s /l P ) 2 



(23) 



so that in order to obtain the elliptical projection from the corre- 
sponding spherical halo we have i) to multiply the overall profile 
by 1/V7 • **) t0 substitute the polar spherical radius with f. The 
intrinsic scale-length has then to be expressed in terms of the pro- 
jected one, see Eq. l |19t . 



4 LENSING INVERSION 

For gravitational lensing studies, the projected map of interest is 
the surface mass density. We will describe the projected NFW den- 
sity in terms of the strength of the lens k b , see Eq. ( 124b . and of the 
projected length scale r s p, i.e. the two parameters directly inferred 
by fitting projected lensing maps. The projected surface mass den- 
sity £ of these density profiles is expressed in terms of the conver- 
gence «, i.e. in units of the critical surface mass density for lensing, 
E cr = (c 2 D s )/(4nGD d D dB ), where D s , D d and D dB are the 
source, the lens and the lens-source angular diameter distances, re- 
spectively. According to our notation in Sec. [3] for a NFW profile, 
the intrinsic l s and the projected lp lengths have to be read as r s 
and r s p, respectively. 

The central convergence of a NFW profile estimated from 
lensing can be written in t erms of C200 and th e projected length 
scale modulus a factor / geo dSereno et al]|2009l) . 



S cr X K s — — ^=/5s^"sP ? 



(24) 



where as usual p s — S c p CI (z) with p CI being the critical den- 
sity of the universe at the cluster redshift. The concentration enters 
Eq. ((24) through 5 C , see Eq. jilt . The estimate of the mass M200 
depends also on the scale-length r s which is known modulus a fac- 
tor \/7/ e A, see Eq. <19t . Then 



M200 = x 200p cr x (c2oor-sp) 3 %|. 



(25) 



In order to estimate M200 and C200 from the projected NFW param- 
eters directly inferred from the lensing analysis, we need to know 
the elongation of the cluster. The problem is intrinsically degen- 
erate and can not be solved based on lensing information alone, 
even in the ideal case of observations without noise. If a cluster is 
elongated along the line of sight, the concentration para meter and 
the virial mass estimate d from lensing are overestimated dGavazzil 
l2005l ; l0gurietalj|2005h . On the other hand, there are more ineffi- 
cient l ensing orientatio ns for a triaxial halo than there are efficient 
ones dCorless et alj|2009l) . 



4.1 Data sample 

We compiled a sample of strong lensing clusters drawing from pre- 
existing lensing analyses. As selection criteria, we retained only 
clusters that are well defined by a single dark matter halo and whose 
lensing data were fitted with an elliptical NFW model. TableQ]lists 
the final cluster sample, together with corresponding NFW param- 
eters and references to where the lensing analyses were performed. 
For clusters that do not have published arc/multiple images red- 
shift, we assumed a source redshift z s = 2.5. Many input data 
were originally presented with asymmetric uncertainties. To ob- 
tain unbiased estimates, we applied correction formulae for the 
m ean and standard d eviation as given by equations (15) and (16) 
in lD'Agostinil d2004l) . Note that there are different definitions for 
the elliptical radius, which affect the numerical value of the pro- 
jected scale-length. We took care to translate published data to the 
notation in the present paper. Furthermore, some studies exploited 
elliptical NFW potential instead of el liptical mass density. When 
necessary, i.e. for the sub-sample from lComerford et alj d2006h . we 
converted the pot ential ellipticity to isod ensity ellipticity according 
to the relation in lGolse & Kneibl d2002l) . As a final precaution, we 
forced errors on re s and r s p to be at least of 10% and the error on e 
to be at least 0.015. Such uncertainties mirror discrepancies among 
different studies of t he sa me data set, see th e analyses of A 1703 in 
Richard et al.l d2009h an d lbguri et alj d2009h or MS2137. 3-2353 in 
Comerford et alj d2006h and lGavazzil d2005h . 

Using the full probability distribution instead of the estimate 
of mean and error for the ellipticity and the central convergence 
would be an improvement. However, we limited our method to 
quite regular clusters (uni-modal and well fitted by a NFW profile). 
From the detailed analyses collected in the literature for each clus- 
ter, we found no evidence for complex parameter distributions, with 
probability functions that are single-peaked and generally well- 
behaved. 



4.2 Inferred parameters 

In order to extract the physical information, i.e. to the determine 
the parameters C200, M200 and eA, we have then to use additional 
constraints. Since we want to test theoretical predictions, we will 
employ the prior from the C200 — M200 relation as given in Eq. (0. 
Such an additional third constraint, together with Eqs. d24l|25t , al- 
lows us to determine the elongation of the cluster and its mass and 
concentration. The prior is very strong so the estimated C200 and 
M200 will fit nicely the theoretical prediction. On the other hand, 
eA is free to take whatever value allows the cluster to fit the lens- 
ing constraints and the C200 — M200 relation at the same time. Un- 
physical values for eA (either < 1 or > 1), that would describe 
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Name 


z d 


Zs 




r sP (kpc/h) 


e 


ref a 




o 

u.zo 


U.OOO 


n 1 q -t- n nzL 


^ZLO -t- QO 


U.J/ u.uoo 




MS 9137 


0.313 


1.501 




112 ± 11 




9 


AC 114 


0.315 


3.347 


22 + 02 


«an + 70 

UOU _1_ 1 VJ 


n ^02 + 01 8 


3 


C1G 2244-02 


0.33 


2.237 


0.18 ± 0.02 


300 ± 30 


0.242 ±0.015 


4 


SDSS J1531+3414 


0.335 


1.096 


1.2 ±0.8 


210 ± 110 


0.47 ±0.23 


1 


SDSS J1446+3032 


0.464 




3.2 ±2.0 


110 ±50 


0.62 ±0.34 


1 


MS 0451.6-0305 


0.55 


0.917 


0.28 ± 0.03 


350 ± 30 


0.425 ±0.015 


4 


3C 220.1 


0.62 


1.49 


0.18 ± 0.02 


320 ± 30 


0.497 ±0.015 


4 


SDSS J21 11-01 15 


0.637 




7.1 ± 1.5 


57 ±11 


0.46 ±0.27 


1 


MS 1137.5+6625 


0.783 




0.26 ± 0.03 


330 ± 30 


0.300 ± 0.015 


4 



Table 1. The st r ong le nsing cluster data sample. References: 1 stands for Ipguri etal] J2009h : 2 for G avazz 1 J2005h : 3 for ISereno et al] 120091) : 4 for 
IComerford et alj £>006). For clusters with multiple image systems, we picked out one source redshift (Col. 2). The central convergence k s for the NFW 
model refers to such redshift; r s p is the projected scale length. 



Name 


C200 


Standard c 2 oo — M200 
A/20o(10 14 M Q /h) 


t-A 




C200 




Enhanced c 2 oo — M 2 oo 
Af 2 oo(10 14 M Q /h) 


eA 


A1703 


2.98 ±0.16 


12 ±4 


0.66 ±0. 


.19 


3.46 ±0. 


.18 


13 ±5 


0.91 ±0.26 


MS2137 


3.41 ±0.13 


2.0 ±0.4 


0.068 ±0. 


Oil 


3.95 ±0. 


.15 


2.3 ±0.5 


0.093 ± 0.014 


AC 114 


2.94 ±0.13 


12 ±3 


1.06 ±0. 


.17 


3.40 ±0. 


.16 


14 ±3 


1.44 ± 0.024 


C1G2244-02 


3.26 ±0.13 


3.3 ±0.7 


0.69 ±0. 


.11 


3.77 ±0. 


.15 


3.7 ±0.8 


0.95 ±0.15 


SDSS1531 


3.1 ±0.4 


7±7 


0.04 ±0. 


.04 


3.6 ±0. 


.4 


8±8 


0.06 ± 0.05 


SDSS 1446 


3.2 ±0.4 


3±3 


0.015 ±0. 


.014 


3.7 ±0. 


.5 


3±3 


0.020 ±0.019 


MS0451 


2.82 ±0.13 


7.6 ± 1.6 


0.29 ±0. 


.05 


3.27 ±0. 


.15 


8.6 ±1.9 


0.40 ± 0.06 


3C220.1 


3.03 ±0.13 


2.6 ±0.6 


0.79 ±0. 


.13 


3.51 ±0. 


.15 


2.9 ±0.6 


1.09 ±0.17 


SDSS2111 


3.0 ±0.2 


2.6 ± 1.6 


0.004 ± 0. 


.002 


3.5 ±0. 


.2 


2.8 ±1.8 


0.006 ± 0.003 


MSI 137 


2.77 ±0.13 


4.3 ±0.9 


0.71 ±0. 


.12 


3.21 ±0. 


.15 


4.9 ±1.0 


0.97 ±0.16 



Table 2. Concentration, mass and elongation for each cluster inferred through lensing inversion assuming as a prior either a standard or an enhanced mass 
concentration relation. Masses are in units of 1O 14 M0 / h. 



more filamentary structures than virialized clusters, will point more 
to outliers with respect to predictions than to extremely elongated 
structures. The most likely explanation for extreme eA values is 
then that the corresponding clusters do not follow the relation im- 
posed a priori. This can be view as a sort of proof ab absurdo. 

Results are listed in Table [2] We cons idered both the C200 — 
M200 as determined in iDuffv et alj d2008h and the case of over- 
concetrated clusters (A = 6.68). To account for measurements er- 
rors, we draw lensing parameters (k s and r s p) from random normal 
distributions with mean and dispersion given by the reported central 
location and scale, see TableQ] Similarly, theoretical uncertainties 
on the mass-concentration relation were accounted for by drawing 
the parameters A, B and C, which describe Eq. {3}, fr om Gaussian 
distrib utions with mean and dispersion values found in lDuffv et alj 
(2008). Then, for each set of parameters (k s , r s p, A, B and C) we 
solved the system of equations Eqs. ll3l l24|[25l >, discarding only so- 
lutions with either c 20 o > 40 or Af 20 o > 1O 18 M //i. The values 
listed in Table [2] are the biweight estimators for l ocation and scale 
of the final inferred distributions jBeers et al . 1990). 

If we use the enhanced C200 — M200 relation, the concentration 
of each cluster increases by ~ 16%, the mass by ~ 13% and eA 
by ~ 27%, i.e. the elongation shrinks. Even if clusters come out 
less elongated if we assume that the lensing population is intrin- 
sically over-concentrated, we see that some outliers are still there. 
Four out of 10 clusters have eA < 0.1, i.e. the size along the line 
of sight should be larger than ten times the maximum length in the 



plane of the sky. Lensing parameters of SDSS 1531, SDSS 1446 
and SDSS 2111 had quite large observational uncertainties which 
propagate in the estimate of the intrinsic cluster parameters. How- 
ever, the estimated values of eA are so small that the ordinary value 
of ~ 1 can be excluded for such clusters at a high confidence 
level. Even doubling the normalization factor of the c — M relation 
(i.e. assuming A ~ 11.4), elongation parameters for two cluster 
(SDSS 21 12 and SDSS 1446) would remain smaller than one tenth 
(e A ~ 0.019 and 0.070, respectively). 

Note that final results on elongation would have been consis- 
tent if we had chosen different methods for deriving the strong lens- 
ing paramet ers. The elongation o f A 1703 calculat ed using the data 
reported in Richard et al. ( 2009) or lLimousin et alj J200S). which 
both fitted the lensing potential, turns out to be 0.3 5 ± 0.11 or 
0.77 ± 0.17 respectively, the value of eA based on lOguri et alj 
(2009), that directly fitted the convergence, being intermediate be- 
tween the two, see Table [2] The elongation of MS 2137 using 
data in Comerfo rd et alj d2006l) , that fitted the lensing potential, is 
0.051 ± 0.008, compatible wit h the result based on direct conver- 
gence fitting in lGavazz il d2005h . see Table [2] Then, independently 
of the lensing technique used, results are quite consistent within 
the errors, both for mildly (A 1703) or very elongated (MS 2137) 
clusters. 

It is quite reassuring that whenever a cluster has been analyzed 
either fitting the potential or the convergence, the estimated elonga- 
tion does not change in a significant way. Together with the central 
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Figure 1. Probability density functions for the elongations of different clus- 
ter populations, see legend. The bars denote the measured values for our full 
sample assuming a standard C200 — relation. 



91,2 


e 


P(e A < 0.1) 


P(e A < 0.2) 


P(e A < 0.3) 


fiat 


flat 


< 10- 3 % 


8 x 10"3% 


2.7 x 10- 3 % 


flat 


biased 


4.5 x 10"3% 


3.5 x 10- 2 % 


1.87 x 10" 2 % 


TV-body 


flat 


< 10-3% 


1 x 10-3% 


9.5 x 10- 3 % 


TV-body 


biased 


< 10-3% 


1.5 x 10"3% 


5.7 x 10- 2 % 



Table 3. Probability (in %) to have an elongation larger than a given thresh- 
old value for different populations of galaxy clusters. 



convergence, our method needs only an estimate of the projected 
ellipt icity, which is quite w ell measured with strong lensing anal- 
yses. iGolse & KneibU2002l) discussed in details how potential and 
surface mass ellipticities are related, and their analysis showed how 
the mass density ellipticity can be estimated using a previous deter- 
mination of the potential ellipticity. Once the ellipticity of the sur- 
face mass density is known, our method relies only on geometrical 
projections and is not affected anymore by lensing non-linearities. 
In fact, we always de-project the surface mass density (instead of 
the potential) to obtain the intrinsic mass distribution. 



5 EXPECTED VS. OBSERVED ELONGATION AND 
ELLIPTICITY 

The chance to observe a very elongated cluster can be assessed on 
a more firm ground. We derived the probability density function 
(PDF) for a given elongation, P(e A ), and a given ellipticity P(e) 
under different assumptions. As discussed in Sec. [3] elongation and 
ellipticity depend on the intrinsic axial ratios, qi and q-2 and the ori- 
entation angles 9 and op. We considered four scenarios. For the axial 
ratios, we considered either the TV-body predictions in Eqs. <[4j [6j 
or a flat distribution, see Eqs. ([7] [8}. For the alignment we con- 
sidered either the biased distribution for P(9) in Eq. ([9} or a flat 
distribution, Eq. dlOb . For the azimuth angle op we always used a 
flat distribution, P(cp) = const. 



5.1 Elongation 

The PDF for the elongation, P(e A ), is plotted in Fig.Q] It is pretty 
evident that populations of clusters preferentially aligned with the 
line of sight make a better job to reproduce the observed sample, 



91,2 


e 


All 


e A > 0.1 


flat 


flat 


4.8 x 10~ 2 


1.5 x 10" 2 


flat 


biased 


2.6 X 10- 2 


1.0 x 10- 2 


TV-body 


flat 


7.9 x 10- 1 


5.9 x 10" 1 


TV-body 


biased 


8.6 x 10"3 


1.2 x 10- 1 



Table 5. Kolmogorov-Smirnov significance level that the measured elliptic- 
ities of the observed samples of clusters (either all of them or the six with 
e A > 0.1) are drawn from a given population. 



apart from the group at e A < 0.1. Cumulative distributions are 
plotted in Fig. [2] and probabilities at very low threshold values are 
listed in Table [3] Chances for very elongated clusters are very tiny. 
Even for biased distributions just one out of few thousands clusters 
is expected to have e A < 0.1. Even in the most favourable case 
of a population of clusters biased in 9 and flat in axial ratios, the 
chance to have four out of ten clusters with e A < 0.1 would be a 
very tiny 1.7 x 10~ 16 . So we can conclude that such clusters are 
very likely outliers of the mass-concentration relations. 

Further quantitative comparisons can be performed exploit- 
ing the Kolmogorov-Smirnov (KS) test. When we consider the full 
sample, none of the investigated populations gives a good fit to the 
data. The better performer, i.e. a population with TV-body like axial 
ratios and biased alignments, give a KS significance level of < 1%, 
both for the standard or the enhanced C200 — TW200 relation, see Ta- 
bleg] 

The significance levels improve very significantly when we 
consider the subsample with e A > 0.1. The prediction from 
TV-body simulations reproduce very well the observed distribu- 
tion both for the standard (~ 31.1%) and the enhanced relation 
(~ 4.7%). For the enhanced relation, also distributions flat in the 
axial ratios perform well for both populations suffering orientation 
bias (~ 2.1%) or unbiased (~ 97.2%). 



5.2 Ellipticity 

The ellipticity distribution of our sample is not near as informative 
as the elongation one. Probability density functions both for unbi- 
ased or biased populations have not negligible values in correspon- 
dence of the observed ellipticities, see Fig. [3] Populations with flat 
axial ratios are preferentially rounder (e > 0, e A ~ 1) since high 
values of q\ are not penalized, but the observed sample do not help 
to discriminate. The KS test is inconclusive too, see Table [5] even 
if the biased TV-body-like population performs remarkably better 
considering the e A > 0.1 subsample. However, the ellipticities of 
such subsample are nothing special. According to a KS test, the 
ellipticities of the outliers (i.e. clusters with e A < 0.1) might be 
drawn from the full sample with a significance level ~ 98%. 

Since our sample is neither homogeneous or statistical, we 
are cautious in drawing conclusions, but some indications seem 
to be quite strong. There is a number of clusters whose over- 
concentration problem can not be solved just considering some par- 
ticular geometrical configurations. Even strong biases in intrinsic 
triaxiality and alignment would not solve the problem. Once such 
outliers are excluded from the analysis, theoretical predictions are 
in very good agreement with data. Populations with an alignment 
bias perform much better than randomly oriented clusters. There 
is also some evidence for intrinsic axial ratios distributed accord- 
ing to the outputs of TV-body simulations, even if, under suitable 
circumstances, flat populations can give good results too. 
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All 


e A > 0.1 


All 


e A > 0.1 


flat 


flat 


9 x icr 8 


2.5 x icr 4 


3 x 10~ 5 


2.1 x 10" 1 


flat 


biased 


9 x icr 6 


3.5 x 10" 3 


9.3 x icr 3 


9.7 x 10- 1 


iV-body 


flat 


6 x icr 7 


2.5 x 10" 4 


5 x icr 5 


9.1 x 10~ 3 


N-boAy 


biased 


7.8 X 10~ 3 


3.1 x icr 1 


9.6 x 10" 3 


4.7 x 10" 2 



Table 4. Kolmogorov-Smirnov significance level that the elongations of the observed samples of clusters (either all of them or the six with e A > 0.1) are 
drawn from a given population. We considered e A obtained from either a standard or an enhanced mass-concentration relation. 
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Enhanced C2oo - ^20i 



u 



m 0.6 



fi 0.4 



i 0.2 



0.0 



0.0 





















flat q, flat 6 






flat q, biased 8 






/V-body q, flat 8 ~ 






■ /V-body q, biased 8 - 



0.5 



1.0 

''A 



1.5 



2.0 



Figure 2. Predicted cumulative distribution functions for elongation versus measurements. The full and dashed step-lines are for the full observed sample and 
for the clusters with e A > 0.1, respectively; the smooth functions are the predicted distributions under different assumptions, see legends. The left and right 
panels show the observed elongations computed assuming either a standard or an enhanced mass-concentration relation, respectively. 



The assumption that lensing clusters are intrinsically over- 
concentrated partially reduce the problem, but expected distribu- 
tions and observations would be compatible with a very low signif- 
icance level of < 1% and the problem of having nearly half of the 
sample with very extreme elongation would be still there. In gen- 
eral, the data analysis performed on our limited sample does not 
provide evidence for intrinsic over-concentrations, an orientation 
bias being enough to account for observations of normal clusters 
(e A > 0.1). 



6 INTRINSIC AXIAL RATIOS AND ORIENTATION 

Knowledge of the sizes of a cluster in the plane of the sky and 
along the li ne of sight a llows us to put constraints on its intrinsic 
geometry JSerendl2007l) . However, even exploiting such strong as- 
sumptions on the shape, inversion can be not unique: intrinsically 
different elli psoids can cast on the plane of the sky in the same way 
dSerendl2007h . In order to infer the properties of the cluster and 
derive its orientation and shape, we have to exploit some external 
information. We will consider two kinds of prior: first a sharp one 
which assumes the cluster to be axially symmetric; then some less 
informative priors on the distribution of intrinsic axial ratios for 
triaxial haloes. 



dPe Filippis et al.l2005l ; lsir"eno et al.l2 006). even if diffuse prolate- 
ness can not be excluded ( Plionis et aljl99ll:lde Theiie et al J 19951: 
Basilakos et alj|2000l : ICooray||2000l : IPlionis et alj|2004 |Paz et all 
2006). Once the elongation of a cluster is known together with its 
project ed ellip ticity, strong constraints can be put on its intrinsic 
shape JSerenoll2007t) . Axial symmetry reduces the number of un- 
known parameters to a couple: a single axial ratio, 1), and the 
inclination angle of the symmetry axis, i. 

A prolate-like solution is admissible when the size along the 
line of sight is larger than the minimum width in the plane of the 
sky, that is, when e A ^ ep. The intrinsic parameters can be written 
as q\ = qi = q and 6 = i. In terms of the measured quantities 
(Sereno 2007), 



e A 



e P 



(26) 



(27) 



An oblate-like solution is admissible when the size along the 
line of sight is smaller than the maximum size in the plane of the 
sky, that is, when e A ^ 1. According to our notation, for an oblate 
ellipsoid, q\ = q, q-z = 1 and cos i = sin 9 sin <p. Inversion gives 
(Sereno 2007), 



6.1 Axial symmetry 

As a working hypothesis, let us first consider if the cluster 
shape can be approximated as an ellipsoid of revolution. Previ- 
ous studies have shown that clusters seems to be quite triaxial 



9 = 



cos I = 



epe A 



(28) 



(29) 



The prolate and the oblate solutions are admissible at the same time 
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Figure 3. Probability function for the projected ellipticity. Left panel: probability density functions for the ellipticity of different cluster populations, see 
legend. The bars denote the measured values for our full sample. Right panel: predicted cumulative distribution functions for a given elongation versus 
observations. The full and dashed step-lines are for the full observed sample and for the clusters with > 0.1, respectively; the smooth functions are the 
predicted distributions under different assumptions, see legend. 







Prolate 






Oblate 




Name 


Comp. 


g 


cos i 


Comp. 


q 


cos i 


A1703 


~ 1 


0.26 ± 0.08 


0.94 ±0.04 


0.0342 


0.59 ±0.05 


0.25 ± 0.11 


MS2137 


~ 1 


0.040 ± 0.007 


0.9994 ± 0.0002 


<10" 5 


NA 


NA 


AC 114 


~ 1 


0.26 ± 0.05 


0.88 ±0.05 


0.622 


0.43 ± 0.04 


0.26 ± 0.09 


C1G2244-02 


~ 1 


0.40 ±0.07 


0.93 ±0.03 


2.31 x 10" 3 


0.74 ±0.03 


0.24 ±0.12 


SDSS1531 


~ 1 


0.012 ± 0.011 


0.999795 ± 0.0002 


<10" 5 


NA 


NA 


SDSS1446 


~ 1 


0.0014 ± 0.0014 


0.999989 ± 0.000015 


<10" 5 


NA 


NA 


MS0451 


~ 1 


0.098 ± 0.017 


0.9902 ± 0.003 


<io~ 5 


NA 


NA 


3C-220.1 


~ 1 


0.20 ±0.03 


0.94 ±0.02 


0.0542 


0.48 ± 0.02 


0.16 ± 0.07 


SDSS2111 


~ 1 


0.0010 ± 0.0009 


0.999998 ± 10" 6 


<10" 5 


NA 


NA 


MSI 137 


~ 1 


0.35 ± 0.06 


0.93 ±0.03 


6.45 x 10~ 3 


0.68 ±0.03 


0.23 ± 0.11 



Table 6. Intrinsic parameters (axial ratio q and inclination angle i) assuming either prolateness or oblateness. The column "Comp." gives the significance level 
for a cluster shape to be compatible with a given set of data. For a very low compatibility with a given shape hypothesis, parameter values are not available 
(NA). 



only when the size along the line of sight is intermediate, i.e. 1 ^ 

eA < e P . 

Results of the inversion are listed in Table [6] Intrinsic param- 
eters have been obtained by means of Eqs. J26M27b for the prolate 
case and Eqs. J28I |29t for the oblate-case. Input values for elon- 
gation and ellipticity were randomly extracted from normal distri- 
butions centred in the measured value and with dispersion equal to 
the observational uncertainty. The values listed in Table [6] are the 
bi-weigth estimators of the final distributions of the inferred param- 
eters. The significance level for a given shape has been obtained as 
the fraction of drawn ep and eA for which a given compatibility 
conditions is fulfilled. We considered only elongation values ob- 
tained assuming a standard mass-concentration relation. The pro- 
late hypothesis is compatible with the full sample but the shapes 
should be extremely long and narrow (q < 0.35) and nearly per- 
fectly aligned with the line of sight (cos i > 0.88). The conclusion 
that clusters with eA < 0.1 are outliers is further stressed by the 
very small space volume allowed for the intrinsic parameters (under 
wrong hypotheses, uncertainties are very likely to be very small). 

A population of oblate clusters do not provide a good descrip- 
tion of the data. Only a tiny region in the parameter space of elon- 



gation and ellipticity allowed by the data is compatible with such 
an hypothesis. Only AC 114, with eA ~ 1, has a good chance 
to be described by an oblate shape (~ 60%), otherwise signif- 
icance levels are < 5%. For the few clusters for which oblate- 
ness is marginally compatible, inclination angles would still be bi- 
ased, symmetry axis being nearly perpendicular to the line of sight 
(cosi < 0.26), whereas intermediate axial ratios would be pre- 
ferred (0.43 < q < 0.74). 

6.2 Triaxial clusters 

In order to exactly determine the intrinsic shape of a triaxial clus- 
ter, we should know both ellipticity and elongation together with 
two additional observational constraints. The problem of invert- 
ing a projected map is intrinsically degenerate and even adding X- 
ray observations or measurements of t he Sunyaev -Zeldovich effect 
would not make the inversion unique dSerencl2007h . An alternative 
appro ach is to use priors on the intrinsic parameters JCorless et al.l 
2009). Here, we try to solve the system of equations 



ca 



e P (qi,q 2 ;9,<fi); 
eA(qi,qr,0,ip) 



(30) 
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Figure 4. Posterior probability density functions for the intrinsic parameters of A 1703. Panels from the left to the right are for the PDF of q\, q2, cos 8 
and cos tp, respectively. Full and dashed thick lines have been obtained assuming a N -body-like and a flat prior on the axis ratios, respectively. The full and 
dashed thin line in the left panel represent the iV-body and the flat prior for P(qi), respectively; the full and dashed thin line in the q^ -panel represent the 
prior distributions according to either a ./V-body or a flat prior, respectively; the thin and dashed full line in the cos 0-panel represent the biased and the flat 
distributions on the orientation angle. Such priors on cos 8 were not used to derive the PDFs. Finally the flat line in the cos ip-panel represents an uniform 
distribution. 
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Figure 5. The same as Fig.|4]for the cluster AC 1 14. 




Figure 6. The same as Fig.|4]for the cluster C1G 2244. 
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Figure 7. The same as Fig.|4]for the cluster 3C 220. 
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Figure 8. The same as Fig.|4]for the cluster MS 1137. 
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TV-body q Flat q 

Name Comp. qi qi cos 9 cos ip Comp. q± q^ cos 9 cos ip 



A1703 


0.094 


0.39±0. 


07 


0.55±0. 


11 


0.87±0. 


.10 


0.81±0. 


16 


0.025 


0.37±0. 


.10 


0.61±0. 


15 


0.90±0. 


09 


0.76±0.17 


MS2137.3 


< 10-4 








NA 








<l -4 








NA 






AC 114 


0.104 


0.36±0. 


06 


0.52±0. 


13 


0.77±0. 


.11 


0.93±0. 


07 


0.054 


0.31±0 


.09 


0.61±0. 


18 


0.84±0. 


12 


0.85±0.13 


C1G2244 


0.149 


0.43±0. 


06 


0.57±0. 


08 


0.93±0. 


.04 


0.67±0. 


20 


0.026 


0.45±0. 


.10 


0.60±0. 


12 


0.92±0. 


05 


0.67±0.22 


SDSS1531 


1.2xl0~ 3 


0.44±0. 


06 


0.60±0. 


10 


0.82±0. 


.25 


0.70±0. 


17 


4.3xl0~ 4 


~ 0.26 




~ 0.49 




~ 0.83 




- 0.60 


SDSS1446 


1.7xl0~ 4 


~ 0.45 




~ 0.70 




~ 0.65 




- 0.39 




1.7xl0~ 4 


- 0.16 




~ 0.76 




~ 0.75 




- 0.63 


MS0451.6 


<10-4 








NA 








3 x 10-4 


- 0.14 




~ 0.27 




~ 0.99 




~ 0.64 


3C220.1 


0.037 


0.33±0. 


06 


0.50±0. 


12 


0.85±0. 


.07 


0.95±0. 


05 


0.025 


0.29±0. 


.09 


0.56±0. 


15 


0.92±0. 


06 


0.86±0.13 


SDSS2111 


<l -4 








NA 








< io- 4 








NA 






MSI 137.5 


0.174 


0.38±0. 


05 


0.52±0. 


08 


0.92±0. 


.04 


0.74±0. 


16 


0.028 


0.41±0. 


.09 


0.59±0. 


12 


0.92±0. 


06 


0.73±0.19 



Table 7. Intrinsic parameters for a triaxial shape (axial ratios, q\ and q2, and orientation angles 9 and tp) inferred using different priors for the intrinsic axial 
ratios distributions. The column "Comp." gives the significance level for a cluster shape to be compatible with a given set of data. 



using a couple of proxies for the axial ratios. In particular we ex- 
ploit either a flat distribution or the guess from TV-body simula- 
tions. Operatively, we randomly extract the axial ratios from the 
assumed prior distribution and then solve for 9 and tp in Eqs. J30| l. 
For each iteration, a couple of input values for ep and eA is also 
randomly drawn. If there is a solution to the system, we consider 
the drawn q± and q2 and the corresponding 8 and tp to be a sample 
from the posterior distribution. 

Results are listed in Table [7] where as usual we reported bi- 
weight estimators. Final estimates are quite insensitive to the priors. 
This asserts the validity of our inversion approach, since a Bayesian 
analysis is effective as far as the effect of priors is as not informa- 
tive as possible. The parameter space for solutions is quite narrow 
for very elongated clusters (eA < 0.3), and we actually were able 
to find very few of them, less than one out of ten thousands draw- 
ings. Assuming the sharp prolate prior on the shape, we could find 
some extremely elongated configurations, but such intrinsic shapes 
are pretty much excluded by assuming more realistic priors on the 
axial distributions as the ones expected for general triaxial con- 
figurations. This further suggests that our sample contains several 
outliers of the mass-concentration relation. 

Posterior probabilities for A 1703, AC 11 4, C1G2244, 3C220 
and MSI 137 are plotted in Figs, g] |5] |6] |7] and [8] respectively. 
The final distributions have been smoothed using a Gaussian ker- 
nel estimator with reflective boundary conditions dVio et alii 19941 ; 
iRvdenll 19961) . For each cluster, whatever the prior on the axial ra- 
tios, the posterior probabilities are quite similar. Even if we assume 
a flat distribution for the axial ratios, the posterior probability is 
quite similar to the prediction from TV-body simulations. Note that 
the alignment bias is confirmed by the above analysis without any 
prior assumption on the orientation. 



7 DISCUSSION 

Recent observational analyses have been finding many lensing clus- 
ters with Einste in radii much larger than expected in a standard 
ACDM model teroadhurst et alj|2008h Ogur i et alj|2009l) . We per- 
formed a statistical analysis on a sample of 10 clusters that were 
well fitted by a single NFW model. Our method was as follow. We 
supposed theoretical expectations from TV-body simulations to be 
true and modelled clusters as NFW haloes fitting standard mass- 
concentration relations. Then, we found the elongation along the 
line of sight of the clusters required to satisfy at the same time 



both lensing data and theoretical predictions. Finally, we studied for 
each cluster which intrinsic shape and orientation were compatible 
with the inferred elongation and the measured projected ellipticity. 
At each step, we checked the exploited hypothesis ab absurdo by 
finding any inconsistency between theoretical predictions and ac- 
tual conditions under which expectations can be in agreement with 
data. 

We first considered the inferred distribution of elongations, 
comparing that with the probability to really see them observ- 
ing a given population. We found two groups in our sample. The 
first group is in very good agreement with what expected from 
a population of clusters fitting the mass-concentration relation 
and preferentially oriented along the line of sight, as suggested 
by theoretical analyses o f lensing clusters dHennawi et al J 120071 : 
lOguri & Blandford 2009). Observed ellipticities and inferred elon- 
gations are also in agreement with intrinsic axial ratios following 
distributions derived in TV-body simulations. There is no evidence 
for such lensing clusters to be intrinsically over-concentrated even 
if data can not exclude that. 

The second subsample in our analysis is made of clusters very 
likely to be outliers of the mass concentration relation. To fulfil the 
expected relation, they should develop along the line of sight as a 
filamentary structure with extreme elongation, a clearly poor de- 
scription for massive haloes. Even allowing for more concentrated 
halos for a given mass by enhancing the c— M relation, elongations 
would still be extreme. 

Even if our sample was not statistical, we took care of se- 
lecting quite regular clusters. However, bimodal structures nearly 
aligned with the line of sight would seemingly have a regular mor- 
phology in the plane of the sky. Such configurations would boost 
the apparent concentration, but they are very rare and it is problem- 
atic to consider all of our outliers within this scenario. Furthermore, 
we used a mass-concentration relation derived for the full sample 
of clusters, not just the virialized and regular ones. 

The second step in our analysis further strengthens such view. 
Using a series of statistical priors, we found for each one of the 
mildly elongated clusters an intrinsic structure and orientation com- 
patible with the inferred elongation. Prolate shapes make a better 
work in explaining data than oblate clusters, but in general data are 
fully compatible with triaxial structures. Whatever the hypothesis 
exploited as prior on the intrinsic axial ratios, inferred intrinsic pa- 
rameters suggests mildly triaxial clusters with an alignment bias, in 
very good agreement with expectations from TV-body simulations. 

An alternative approach would have been to apply the fitting 
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procedure to inclined triaxial haloes in the first place. As far as the 
estimates of the projected ellipticity and of the central convergence 
are not biased, our method should be able to fully explore the space 
of the triaxial parameters (shape and orientation). In fact, given a 
projected map, we consider all the intrinsic configurations com- 
patible with data. This is done through Eqs. d 141 124t . Such set of 
equations allows us to study the full intrinsic parameter space, even 
in case of disjoint regions compatible with data. 

Either projecting intrinsic parameters and fitting to the lensing 
data or fitting projected maps and then deprojecting (as done in this 
paper) should give the same final result. In fact, in a pure lensing 
analysis the triaxial structure of a cluster is constrained only by its 
projected map, so that both procedures should pick out the same 
sets of intrinsic parameters that can fit the measured quantities, i.e. 
ellipticity, orientation and central surface density. 
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